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@ Method and system for classifying agglutination reactions. 

(g) A method and system for analyzing a solutton for an agglutination patlem. The method comprises the 
steps of producing an illuminated image of the solution on an an^y of pixels, and assigning to each 
pixel In the Illuminated image, a data value representing the intensity of the Siumlnated Image on the 
pixel. With the preferred processing procedure, the pixel anray Is separated Into a plurality of zones, and 
the data values for the pixels in each zone are processed according to a respective predetennlned 
procedure to determine values for a predefined set of variables. Then, those determined values are 
processed to determine whether an agglutination pattern is present in the solution, and if so, to classify 
that pattern Into one of the predefined classes. 
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BACKGROUND OF THE INVENTION 

This Invention generaily relates to a method and systemfor the detection and quantification of aggiutinates 
formed in response to immunological agglutination reactions, and more particularly, to such a method and sys- 
5 tern employing automatad image and data procesaing means to automatically detect and classify agglutination 
patterns. 

immunological agfllutinatlon reactions are used for Identifying various kinds of blood types and for detect- 
ing various kinds of antibodies and antigens In blood samples and other aqueous solutions. In a conventional 
procedure, a sample of red blood cells is mixed with serum or plasma in test tubes or mlcroplates, and the 

10 mbcture may then be incubated and centrif uged. Various reactions either occur or do not occur depending on, 
for example, the blood type of the red blood cells or whether certain antibodies are present in the blood sample. 
Typically, these reactions manifest themselves aa clumps of cells or particles with antigens and antibodies on 
their surfaces, raferred to as agglutinates. Thus, the absence of any such clumps indicates that no reaction 
has occurred; and the presence of such clumps indicates that a reaction has occun^d. with the size and amount 

15 of such clumps being a <|uantltative Indicator of the level or concentration in the sample, or an indicator of the 
reaction strength, affinity of the complex for which the blood sample was tested. 

Recently, a new agglutination test method -refen-ed to as coiumn agglutination technology, or CAT- has 
been developed. Coiumn Agglutination Technology may be defined as the analysis of blood and blood products 
utilizing filtration as a means of separating agglutinated, precipitated, absorbed, or adsorbed particulate conv 

20 ponenta from non-reactive oomponenta for innmunoassay applications. In this method, gel or glass bead ml- 
croparticles are contained within a small column, referred to as a microcolumn. A reagent such as antl-lgG is 
dispensed in a diluent in the microcolumn and test red blood ceils are placed in a reaction chamber above the 
coiumn. The column, which is typically one of a multitude of columns formed in a transparent cassette. Is cen- 
trif uged. The centrifuging accelerates the reaction, If any, between the reagent and the blood cells, and also 

25 urges the cells toward the bottom of the column. The glass beads or gel in the microcolumn act as a filter, how- 
ever, and resist or Impede downward movement of the particles in the column. As a result, the nature and dis- 
tribution of the particles In the. microcolumn after centrifuging provides a visual Indication of whether any ag- 
glutination reaction occunred In the microcolumn, and If so, of the strength of that reaction. 

In particular, if no agglutination reaction occurs, then ail or virtually ail of the red blood ceiis in the micro- 

30 coiumn pass downward, during centrifuging, to the bottom of the column and form a pellet at that bottom. If 
there Is a very strong reaction between the reagent and the red blood cells, vIrtuaDy all of the red blood cells 
agglutinate, and large agglutinates form at the top of the microcolumn, above the gel or glass beads contained 
therein. The gel or glass beads prevent the agglutinates from passing, during centrifuging. to the bottom of 
the column, so that after centrifuging the agglutinates remain above the gel or beads. 

35 If there Is a reaction, between the reagent and the blood cells, but this reaction is not as strong as the above- 
descri bed very strong reaction, then some but not all of the red blood ceils agglutinate. The percentage of red 
blood cells that agglutinate and the size of the agglutinated particles both vary directly with the strength of 
the reaction. During centrifuging, the unreacted blood cells pass to the bottom of the column, and the distance 
that the agglutinated particles pass downward through the column depends on the size and number of those 

40 particles. Hence, the size of the pellet of red blood cells at the bottom of the microcolumn, and the extent to 
which the aggiutinates penetrate Into the get or glass beads In the mterocoiumn, are both inversely related to 
the strength of the reactton between the reagent and the red blood cells. 

With this CAT, after the desired processing steps have been perfonmed, the microcolumn is observed, or 
read, by a human operator, who then classifies the reaction between the reagent and the red blood cells. Con- 

45 ventlonally, the reaction Is classified as either negative or positive; and If positive, the reaction Is then further 
classified into one of four classes depending on the strength of the reaction. A highly skilled operator is needed 
to properly read and classify the reaction. 

SUMMARY OF THE INVENTION 

60 

An object of this Invention is to automatically analyze aqueous solutions for agglutination patterns. 

Another object of the present invention is to automatically read and classify agglutinatton reactions be- 
tween red blood cell antigens and antibodies. 

Af urther object of this Invention is to produce an image of a blood sample and to analyze that image using 
66 high speed Image and data processing equipment to determine If the blood sample contains an agglutination 
pattern and, if so, to classify that pattern. 

A still another object of the present invention is to provide a system and method for automatically reading 
and classifying agglutination reactions that occur in a column having a microf liter that produces different ag- 



2 




EP 0637 744 A1 

glutination patterns depending on the strength of that reaction. 

Another object of this invention is to provide a system for reading and classifying blood samples thai may 
be used alone or as part of an Integrated, fully automated blood bank system. 

Afurther object of this invention is to provide an accurate, high speed method and system for automatically 

5 reading and classifying red blood cell agglutination reactions. 

These and other objectives are attained with a method and system for analyzing a solution for an agglu- 
tination pattern. The method comprises the steps of producing an illuminated image of the solution on an array 
of pixels, and assigning to each pixel in the Illuminated image, a data value representing the intensity of the 
Illuminated image on the pixel. Those data values are then processed according to a predetermined program 

10 to determine If an agglutination pattern Is present and, if so, to classify that pattern Into one of a plurality of 
predefined classes. With the preferred processing procedure, the pixel anray is separated into a plurality of 
zones, and the data values for the pbcels in each zone are processed according to a respective predetenrnined 
procedure to detemnlne values for a predefined set of variables. Then, those determined values are processed 
to determined whether an agglutination pattern is present in the solution, and if so, to classify that pattern 

f 5 Into one of the predefined classes. 

With the preferred embodiment of the invention disclosed herein in detail, the solutions are contained in 
a column having glass microbeads. The Image processing program searches the location of the column In the 
source image on the pixel anray; and after the column is located, the program generates a window to cover 
the column where the red ceils are located. The program then selects three reference regions from inside and 

20 outside the column and measures the intensity or gray levels In these regions, and these measured gray levels 
are used to determine certain threshold values that are subsequently used In the processing program. 

The cell pellet ia extracted by applying global threshold values in a V-shaped, lower portion of the column, 
and parameters related to the shape of the cell pellet are also calculated. The program then generates a fbced 
mask to cover the bead area in the column. 

25 For the feature calculation, the bead column is divided into five different zones. The region on top of the 
bead column Is defined as the positive zone, the region at the bottom of the column Is defined as the negative 
zone, and the area between the positive and negative zones Is separated Into three intermediate zones. The 
red cells located In the positive zones are extracted using a Uireshold method, and the red cell aggluUnates 
located in the intermediate zones are extracted using a morphological filter. In addition, the balance of the red 

30 cells between the left and right sides of the column Is determined. For each column, the above parameters 
are preferably calculated for both front and back side images of the column, and the two calculated values for 
each parameter are combined. The agglutination reaction ia then classified on the basis of these combined 
features. 

Further benefits and advantages of the Invention will become apparent from a consideration of the follow- 
3S ing detailed description given with reference to the accompanying drawings, which specify and show preferred 
embodiments of the invention. 

BRIEF DESCRiPTION OF THE DRAWINQS 

40 Figure 1 is a schematic diagram of an automated blood analysis system embodying the present invention. 

Figure 2 la a block diagram illustrating several components of the system of Figure 1 . 

Figure 3 Is a front view of a solution cassette that may be used In the system of Figure 1 . 

Figure 4 is a aide view of the cassette. 

Figure 5 Is a top view of the cassette. 
45 Figure OA shows a glass pattern that may be used to focus the camera of the system shown In Figure 1 . 

Figures 6B and 6C show two pattern signals that may be produced on the camera using the glass pattern 
of Rgure 6A, depending on whether the camera Is in focus or out of focus. 

Figure 7 Is a more detailed drawing of the processing subsystem of the analysis system of Figure 1 . 

Figure 8 schematically illustrates the memory board of the image processes of the processing subsystem. 
SO Figure 9 la a front view of the transport subsystem of the analysis system shown in Figure 1 . 

Figure 10 Illustrates the transport subsystem In the analysis system. 

Figures 11 A-11 E show different agglutlnatton patterns that nmy be produced In a column of the cassette 
shown In Figures 3-5. 

Figure 1 2 generally outlines the preferred procedure for processing the Image data produced In the system 
66 of Figure 1 . 

Figure 1 3 Illustrates one step in identifying the locations of the image of the columns in the pbcel array of 
Figure 1. 

Figure 14 lliuatrates a step in identifying the edges of a column Image. 
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Figure 15 shows various reference areas on the pixel array that are used to determine a set of reference 
values. 

Figure 1 6 shows a pellet of red blood cells at the bottom of a column. 
Figure 17 shows a mask used In the image processing. 
5 Figure 18 illustrates the column separated into multiple zones. 

Figure 19 Illustrates the results of a top-hat transformation to a line slice of an Image. 
Figure 20 shows the two parts of the column used to determine the balance of the red blood cell agglu- 
tinates In the column. 

Figure 21 is a different version of the Decision Tree. 

10 

DETAILED DESCRIPTION OF THE PREFERRED EMBODIMENTS 

Figures 1 and 2 Illustrate automated optical reading system 10. generally, comprising holding means 12. 
illumination means 14, Imaging subsystem 16, and processing subsystem 20; and preferably system 10 further 

19 Includes transport subsystem 22, storage means 24, waste receptacle 26, and bar code reader 30. With the 
embodiment of system 10 shown In Figure 1, holding means 12 Includes base 32 and frame 34; and Illumination 
means Includes a pair of fluorescence lights 36a and 36b, diffuser 38, and neutral density filter 40. Imaging 
subsystem 16 includes pixel array 42, housing 44, and lens assembly 46; and this lens assembly, in turn, in- 
cludes lens 50, filter 52, and lens housing 54. Also, the preferred processing subsystem 20 includes prepro- 

20 cessor S6, main processor 60, and Input means such as keyboard 62: and the preferred transport subsystem 
22 shown in Figure 1 1ncludes support means 64 and mover 66. 

Generally, holding means 12 is provided to hold a test sample for analysis, and illumination means 14 is 
provided to produce an illuminated image of the test sample on imaging subsystem 16. Subsystem 16 gener- 
ates a set of signals representing the Illuminated Image formed thereon and then transmits those signals to 

26 processing subsystem 20. The processing subsystem receh/es those signals from subsystem 1 6 and process- 
es those signals according to a predetermined program to determine whether an agglutination pattern is pres* 
ent In a test sample being analyzed and, if so, to classify that pattern In one of a plurality of predefined classes. 

The preferred embodiment of system 10 described herein In detail is particularly well suited for analyzing 
blood samples, and these samples are often referred to as solutions. It should be noted that the present in- 

30 vention may be embodied in systems that analyze other materials, Including other aqueous solutions such as 
urine. It Is not necessary, though, that the material being analyzed be a liquid or a fluid; and. thus, the term 
"solution" as used herein Is used in the general sense as any mixture of liquid, gaseous, or solid substances. 

In addition, the test samples analyzed In system 10 are preferably held within containers, and a large va- 
riety of types and sizes of containers may be used with system 1 0. However, the preferred embodiment of sys- 

35 tem 10 described herein In detail is particularly well suited for use with containers of the type shown at 80 in 
Figures 3, 4, and 5. These containers, referred to as cassettes, are made from a transparent, integrally molded 
plastic material. A multitude of cavities or wells 82, refenred to as columns or mlaocolumns, are formed in the 
cassettes and extend downward from the top edge of the cassette, and, for example, the cassette shown in 
Figures 3-5 contains sbc such microcolumns. 

40 A multitude of very small, transparent glass beads, having diameters on the onjer of magnitude of 1 0 to 
1 00 micrometers, are deposited in and fonm a filter in the lower portion of each microcolumn. Alternately, the 
lower portion of each microcolumn may be provided with a suitable gel that functions In the same general way 
as the mterobeads. Reagents may be pre-dlspensed In the columns of the cassette; and after the columns of 
the cassette are provided with the desired materials, a foil 84 Is typical ly secured on the top edge of the cassette 

4S to cover and close the tops of columns 62. 

When any particular cassette 80 is used, one, some, or all of the mio^ocolumns 62 in the cassette may be 
used. Moreover, each cassette may be used with blood samples from one or more Individuals. In each micro- 
column that is used, a sample of red blood ceils and one or more reagents that react with known agents are 
pipetted into the microcolumn to test that blood sample for the presence of that one or moTB agents. The cas- 

50 sette may be incubated and is then centrif uged. If an agent, for which the blood sample is being tested, is pres- 
ent In the mlcrooolumn. the agent reacts with the red blood cells to fonm agglutinates; and the number, size, 
and distribution of the agglutinates In the microcolumn Is an Indication of the strength of that reaction. 

With reference again to Figure 1 , frame 34 of holding nneans 1 2 forms an elongated channel 34a for holding 
a test sample such as cassette 80; and preferably, as shown In Figure 1 . the longitudinal ends of the channel 

66 34a are open to facilitate or to allow sliding movement of the test sample Into, through, and then from channel 
34a. Also, frame 34 is preferably rotatably mounted on base 32 for pivotal or rotary movement about a central 
vertical axis, and a motor Is connected to frame 34 to pivot or rotate the frame about that axis. 

Illumination means 14, which preferably comprises a pair of fluorescence lamps 36a and 36b. directs light 
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through the test aample held in frame 34 end onto imaging fiubey&tem 16, end epedf icelly onto pixel errey 42, 
which then generates e series of signals representing the test sample. More particularly, pixel array 42 is dis- 
posed inside a camera housing 44, and the pixel array Is preferably comprised of a multitude of light sensors 
each of which is capable of generating a respective one electric current having a magnitude proportional to 

s or representing the intensity of light incident on that sensor. Preferably, these light sensors, or pbceis. are ar- 
ranged in a uniform grid of a given number of uniformly spaced rows and columns. 

With reference again to Figure 2, lens 50 and filter 52 are located forward of pixel anray 42 and are coaxlally 
aligned with each other and with the pixel array, and lens 50 is positioned so that the pixel anray is at the bade 
focal plane of this lens. Preferably, lens 50 and filter 52 are mounted Inside housing 54, which in turn Is nrMunted 

10 on the front end of camera 44. 

As will be understood by those of ordinary skin in the art, any suitable light source 14, lens 50, filter 52, 
and camera 44 may be used In system 10. For examiple. In an embodiment of system 10 that has been actuaDy 
reduced to practice, canrwra 44 Is a Sony XC-75CE video camera, and the pixel anray or sensoring element In 
this camera Is a charged coupled device (CCD) comprising a matrbc of pbceis in a rectangular anray. 752 pixels 

15 by 582 pbceis. The distance between the camera and the cassette held In frame 34 was adjusted so that each 
image on the pixel array contains two columns 82 of the cassette, and the width of each column In the image 
is about 152 pbceis. 

AComponon microiens manufactured by Schneider Corporation was set at an F stop of F/4.0 and attached 
on the camera via an adaptor. Between the lens and the CCD element was f bced a band pass fitter with a center 

20 wavelength of 550 nm and a band vridth of 40 nm. This f llterenhances ths image of the red cells and Improves 
the signal to noise ratio, and the filter was selected on the basis of a spectro-photonneter measurement, which 
Indicates that red blood ceils have Increased absorption of light In the corresponding wavelength range. 

The camera 44 was focused using a piece of glass printed with special patterns. These patterns are dark 
vertical lines with different sizes and widths as shown In Figure 6A. if the camera is in focus, the signal profile 

25 of the pattern (mage is a square wave, as shown in Figure SB. If the camera is out of focus, the signal profile 
looses the sharp edge, as shown in Figure 6C. The pattern signal Is derived, and the optimal focus is reached 
when the maximum derivative value is obtained. 

In this embodiment of system 10 that has been actually reduced to practice, light source 14 Indudes two 
constant fluorescence tubes 36a and 36b (manufactured by Phillips, PL-S, 7-W), a diffuser 38, and a neutral 

30 density filter 40. As partteularly shown In Figure 2, one fluorescence tube 36a was mounted in the front of the 
cassette, and the other fluorescence tube was located upward and rearward of the cassette. The neutral den- 
sity filter 40 is located below the fluorescence tube 36b, and this filter is used to reduce the amount of light 
transmitted to the cassette from tube 36b. The diffuser 38 Is located parallel to and approximately 1.0 mm 
forward of cassette. 

35 Figure 7 is a block diagram Illustrating processing subsystem 20 in greater detail. In this subsystem, the 
electric signals from the pixel anray in camera 44 are conductecJ to preprocessor 56, which may be, for example, 
an image processing board made by Euresys S.A. of Belgiunn. This image processor then converts the electric 
signal from each pbcel of array 42 Into a respective one digital data value and stores that data value at a memory 
location having an address assodated with the address of the pixel that generated the electric signal. The elec- 

40 trie signals being transmitted to Image processor 56 may be Identified in any suitable way with the specific 
pixel that generated the signal. For Instance, the signals from the pbceis of anray 42 may be transmitted to the 
image processor In a given, timed sequence, and a clock signal may be transmitted to the Image processor 
from the camera to Identify the start, or selected intervals, of that sequence. Alternately, each signal trans- 
mitted to the Image processor may be provided with a header or another data tag identifying the particular 

45 pbcel that generated the signal. 

The Euresys image processing board consists of 1 MB memory. As shown in Figure 8. the memory is div- 
ided into four quads; QAU, QAD, QBU, and QBD. Each quad contains a frame of 512x512 pixels. The source 
Image data, S, Is located In one quad QAD, and Quad QAU contains two fbced masks used for covering the 
two columns in the later image processing. Quad QBD is used for a procedure referred to as top-hat processing, 

80 discussed below, and Quad QBD is used fbr temporary operation. 

The data values stored in Image processor 56 are available to main processor 60, which Is connected to 
the Image processor to obtain data values from and to transmit data values to that Image processor. As ex- 
plained in greater detail below, processor 60 is programmed to process and to analyze the data values stored 
in the image processor to Identify the agglutination pattern, if any. in the test sample being analyzed. 

66 Preferably, main processor is, or Is a component of, a personal computer also having keyboard 62 and 

terminal 64. Keyboard 62 is connected to processor 60 to allow operator input thereto, and tenninal 64 is used 
to display visually data or messages being Input Into the processor. In addition, monitor 66 may be connected 
to processor 56 to produce video Images from the data value stored In the processor or In Image processor 
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56. For example, the 6 data values may be tranemitted to monitor 66 to produce thereon an image of the real 
Image produced on pixel array 42. Other eets of data values may be transmitted to the monitor to produce 
refined or processes Images of the real Image. Printer 68 may be connected to processor 60 to provide a visual, 
permanent record of selected data values transmitted to the printer from the processor. 

5 As will be understood by those of ordinary skill in the art, subsystem 20 may be provided with other or 

additional Input or output devices to allow an operator or analyst to interact with processors 66 and 60, Also, 
the Individual components of subsystem 20 are conventtonal and well-known by those of ordinary skill in the 
art. For example, nwnltor 66 may be a high resolution color monitor, and, as discussed above, processor 60 
may be a personal computer, and image processor 56 may be an image processing board made by Euresys 

10 SJ^ of Belgium. 

With reference again to Figure 1 , storage means 24 is located adjacent holding means 12 and is provided 
for holding a multitude or test samples, and preferably indexing means such as a stepper motor is provided 
for moving the storage means through a series of positions to align each of the test samples held therein with 
the holding means. The storage means 24 shown In Figure 1 is particularly designed for holding cassettes 80, 
IS and the storage means fonms a multitude of channels or slots 24a for holding those cassettes. The Indexing 
means moves this storage means 24 so as to al^n each of the channels 24a with the elongated channel 34a 
of frame 34, allowing the cassettes to be slid from the storage means and Into that frame. 

With the preferred embodiment of system 1 0 shown in Figure 1 , storage means 24 comprises a rotatable 
carousel including a rotatable base and a multitude of compartments. Each compartment forms a respective 
20 one of the channels or slots 24a. and each of these slots extends along a radius of the carousel. Further, the 
Indexing means may comprise a stepper motor, and each time the motor is actuated, the motor moves the car- 
ousel so as to align one of the slots 24a with channel 34a of frame 34. Thie stepper motor may be operated, 
for example, to rotate carousel to align one slot 24a at a time with channel 34a, in a clockwise or counterclock- 
wise sequence around the caroussl. Alternatively, the stspper motor may be provided with a programmed, or 
26 programmable, controller that operates the stepper motor to align slots 24a with channel 34a according to that 
program and In an order that may depend on a multitude of variable factors. 

Waste receptacle 26 Is provided for receiving the test samples from holding means 12 after the desired 
imaging has been completed. For example, the waste receptacle may be a container located below and adja- 
cent the output end of channel 34a of frame 34, and positioned so that the test samples that are slid out from 
30 channel 34a fall into roceptacle 26 under the force of gravity. 

Transport subsystem 22 is preferably provided to move test samples, particularly cassettes 60. into and 
then from holding means 12, specifically frame channel 34a. More particularly, with reference to Figures 1. 
9, and 10, support means 64 supports nrtover 66 for sliding movement between carousel 24 and waste recep- 
tacle 26 and over frame 34. In use, mover 66 Is positioned over the carousel, and as the carousel rotates to 
55 align a cassette with frame slot 34a, that cassette is movsd into engagement with the mover. The mover then 
slides the cassette from carousel 24. into frame 34 and into a position directly forward of pbcel array 42. After 
the desired imaging of the test sample is compieted, the mover 66 is operated to slide the test sample through 
the output end of channel 34a and Into wests receptacle 26. Alternatively, depending on the results of the ana- 
lysis of the test sample, that test sample may be moved back into carousel 24, or to another location where 
40 the test sample may be stored, for example, for further tests or for analysis by an operator. 

With particular reference fa> Rgures 9 and 1 0. support means 64 includes a horizontal bar 92; and this bar 
extends from a position directly over carousel 24 to a position directly over waste receptacle 26. and bar 92 
is supported in any suitable manner. Mover 66, In turn, is supported by and is mounted on bar 92 for sliding 
movement therealong, and a suitable motor or motors (not shown) are provided to operate the mover. Mover 
4S 86 may be operated in response to electric signals received from sensors or timers or both, to move the test 
sample in the desired manner. Aiternatlveiy, the mover may be controlled by a programmed, or programmable, 
processor that operates transport subsystem 22 in a predetermined manner and according to a multitude of 
variable factors. 

Preferably, each cassette 80 is provkJed with a bar code 86 identifying selected data about the cassette, 
80 and bar code reader 30 is provided to read the bar code on each cassette and to transmit the data thereon to 
processor 60. For example, the bar code on the cassette may IdenUfy the cassette typo, the date of manu- 
facture of the cassette, and a recommended expiration date for the cassette. The bar code may Include other 
data that identify the cassette manufacturer as well as the time and place of manufacture. As shown in Figure 
1 , the code reader, which may be a standard bar code reader, Is preferably located between storage rack 24 
65 end frame 34 so that the reader scans the bar code on each cassette as the cassette is being transferred from 
the storage rack and into thef rame 34. As an option, If the barcode 86 does not properly idenOfy all the selected 
data, system 1 0 may be operated so that there Is no processing of any Image data from the cassette 80. For 
example, this may be done by not producing any image of the cassetta on pbcel array 42, or If an Image Is 
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produced, by not processing that image. 

In the operation of system 10, a multitude of test samples are piaced in carousel 24. and the carousel Is 
rotated to align a selected one of the slots 24a with channel 34a. Then, nnover 66 slides the test sample In that 
selected carousel slot, Into the desired position In frame 34, and illumination means 14 then directs a beam 
5 of light through the test sample and onto pixel anray 42. Cassette 80 may include positioning marks 88 between 
the columns 62. or elsewhere on the cassette, to help align the cassette In frame 34; and system 10 may be 
operated so that. If the cassette does not have any such positioning marlcs, or If the marks are not properly 
aligned In frame 34, the cassette is rejected or there Is no anaiysis or processing of the cassette. 

Each pbcel of array 42 generates a respective one electric cun-ent having a magnitude representing the 
10 intensity of the light incident on that pixel; and these output currents are converted to digital data values, re- 
ferred to as gray values or as the S values, and stored In image processor 56. Preferably, frame 34 is then 
rotated ISO*", and illumination means directs another light beam through the test sample to produce a further 
Image of the test saniv>Ie on the pixel array. Each pbcel of array 42 generates another respective one electric 
current having a magnitude representing the Intensity of the light of the second Image that Is Incident on the 
15 pixel. These output currents are converted to a second set of digital data values, which are also stored in the 
image processor. These two images of the test sample that are produced on the pbcel array are referred to as 
the front and back images, respectively. 

Processing subsystem 20 then analyzes the images produced on the pixel array, In a manner discussed 
in detail below, to classify the particle patterns in the test sample, and on the basis of that classification, to 
20 classify the reaction between the reagent and the blood cells In the test sample. After the Image processing 
Is complete, mover 66 may be used to slide the test sample Into waste receptacle 26. Alternatively. If processor 
60 detemnines that the test sample should receive further analysis or, for some reason, should be specifically 
brought to the attention of a human operator, then the test sample may be carried, by mover 66, another mech- 
anism, or an operator, to a separate holding area. 
26 in the meantime, carousel 24 is rotated to align a second carousel slot with frame channel 34a; and once 
mover 66 is available, that mover may be used to move the test sample from that second carousel slot into 
the frame channel. The desired Image or Images of that second tost sample are produced and then processed 
to classify the reaction between the blood cells and the reagents In the test sample. 

The above procedure may be continuously repeated until stopped by an operator, or an automated con- 
so troller may be provided to stop the procedure in accordance with a predetermined program. During the oper- 
atton of system 1 0, new test samples may be placed in carousel 24 either by an operator or by a suitable mech- 
anism. 

System 10 is particularly well suited for use in an automated solution testing system or instrument For 
example, a blood anaiysis system or Instrument In which system 10 may be used, is disclosed in contempora- 
35 neously filed patent applicatfon No, , for "An Automated Blood Analysis System", the dis- 

closure of which is herein incorporated by reference (claiming priority frem USSN 075303, Attorne/s 
Ref./P12276EP). 

As previously mentioned, when the test sample Is held In a column 82 of a cassette 80, the number, size, 
and distribution of particles in the column is an indication of whether an agglutination reaction occurred in that 

40 column and, if so, of the strength of the reaction. Conventionally, the reaction is classified as negative (if no 
reaction occunred) or as positive (If a reaction has occurred), and if positive, the reaction is further classified 
as a dass •I'l, *2, *3, or 44 reaction depending on the strength of the raactton. 

Figures 11A through 11E Illustrate these fh^e types of reactions -negative, dass +1. dass +2. dass +3, 
and dass +4, respectively- In a CAT using cassettes 80 of the type shown in Figures 3 and 4. With reference 

45 to Figure 11 A. In the case of a negative reaction, the red blood cells do not agglutinate and, during centrlf uglng, 
those cells pass to the bottom of the column and fonma pellet 102 at that bottom. In a weak, or dass +1 , positive 
reaction, some of the blood cells agglutinate and fomn a relatively small number of small agglutinated partides 
1 04; however, most of the red blood cells do not react During centrlf uglng, agglutinated particles become dis- 
tributed In the lower half of the microbead column, and the unreacted red blood cells pass to the bottom of the 

50 column and fbnm a pellet 1 06 thereat that is slightly smaller than the pellet 1 02 formed in the case of a negative 
reaction. 

In a dass *2 reaction, which Is slightly stronger than a class +1 reaction, a larger percentage of the red 
blood cells agglutinate and the agglutinated partides that fonm are largen however, an appredable number of 
the red blood cells still do not react With reference to Figure 11C, during centrifuging, the unreacted blood 
55 cells pass through the column bottom and form a small pellet 110, and larger agglutinate partides 112 become 
distributed throughout the length of the column of microbeads. in a dass +3 reaction, most or virtually all of 
the red blood cells agglutinate and the agglutinated partides that form tend to be larger than those that form 
in a dass *2 reaction. As shown In Figure 11 D, even after centrifuging, most of the agglutinatad partides 114 
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remain in the upper half of the bead column. In a dass +4 reaction, all, or virtually all. of the red blood cells 
agglutinate and form large agglutinates 116 at the top of the glass beads. These beads prevent the agglutinates 
from passing downward during centrifuging so that after centrifuging, the agglutinates remain above the 
beads, as shown in Figure 11E. 

6 The reagent and the glass beads In column 82. as well as the body of cassette 80. ara substantially trans- 
parent; however, the agglutinated particles and the red blood cells are only partially transparent Hence. In the 
operation of system 1 0. when light Is transmitted through column 82 and onto pixel anray 42, the portion of the 
light that passes through the agglutinated particles and the red blood cells Is incident on the anray at a relatively 
low intensity, while the rest of the light passing through column 82 is Incident on the array at a higher intensity. 

10 Accordingly, in the Image of column 82 that Is formed on pixel array 42, the agglutinated particles and the red 
blood calls appear as gray or shadow amas relative to the rest of the Image. 

With reference to Figure 12, the preferred procedure for processing the image produced on pixel array in- 
cludes four programa- (1) Image acquisition, (2) column detection, (3) feature extraction, and (4) reaction clas- 
sification. The Image acquisition program relates to the management of and the Interface between camera 44 

15 and preprocessor 56, and the column detection program identifies the borders of the columns that appear on 
the pixel array. The feature extraction program extracts Information related to the agglutination reactions from 
the source Image and translates that Information Into quantitative data. The reaction classification program 
separates the reactions into different classes on the basis of the extracted features. 

The first step in the Image acquisition progmm is an initialization step, in this step, the video memory is 

20 Cleared, several variables, discussed below, are set to associated values, and camera 44 is placed In checked 
mode. Then, after a cassette 80 is positioned In front of the camera by the transport subsystem 22. the Image 
acquisition program sends a synchronization signal to the camera, and the camera shutter operates so that 
an Image of the cassette is produced on pixel array 42. The Image data values are then obtained by the image 
processing platform and converted Into digitized signals that are stored in the Image. processor. To optimize 

26 the source image signals, the gain and offset of the electronic board is preferably adjusted with a gray scale 
so that the gray level of a black strip is zero, and the gray level of the glass beads is 170. 

After the Image acquisition program Is completed, the column detection program begins. The first step In 
this program Is to search for the positions of the two columns In the image frame, and in particular, to search 
for the left, right and bottom edges of the two columns on the pixel array. This Is done by using the fact that, 

30 due to the diffraction of light by those edges, those edges appear on the pbcel array slightly darker than the 
Immediately acUacent areas. More specifically, with reference to Figure 13, the program creates two rectan- 
gular areas on the pbcel array to cover the regions where the two columns are expected to be located. The 
size of each araa may be, for example, 220 by 500 pbcels, and the position of each rectangle Is determined by 
preset x and y coordinates of the upper left comer of the rectangle. These parameters are listed below: 

35 



40 



45 



Name 


Value 


Description 


COL LOCA WIDTH 


220 


width of the frame 


COU LOCA HEIGHT 


500 


height of the fr/in» 




1 COUI LOCA <3RGX 


i ^ 


x-coorainate of tlie left column 


C0L2 i,OCA-OUGX 


1 292 


x-c;oordina»:e of the right column i| 



The Y-coordinates for the upper left corners of both areas are zero. 
» To detect the left and right edges of the columns, the program then generates three small areas within 
each rectangular window, as shown In Figure 14. The positions of these areas are defined by the following 
coordinates: 



65 
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Name 


Value 


Description 


UPPER_BORDER_LENQTH 


25 


segment length 


UPPEF^BORDER^THICKNESS 


199 


projeotion thickness 


UPPER_BORDER_Y 


160 


y-coordinate of the middle of the area 


LEFT.BORDER^X 


30 


x-coordlnate of the middle of the left area 


RIGHT,BORDER.X 


190 


x-coordlnate of the middle of the right area 


BOTTOM.B0RDER_LENGTH 


45 


segment length 


BOTTOM_BORDEI^THICKNESS 


19 


projection thickness 


BOTTOM_BORDEFLY 


435 


y-coordlnate of the middle of the area 



The two left and right symmetric areas are used to detect the left and right side borders of the column. 
The gray values In the blocks are projected into two one-dimension vectors. At each point on each vector, a 
projected value Is obtained by adding the gray values at all pixels, in the corresponding vertical line through 

20 that point The length of the vector is defined by the UPPER_BORDER_LENGTH. These values are then de- 
rived and the maximal derivative is obtained. This maximal derivative conresponds to the maximal variation of 
the gray scale, and on the basis of Its location, the edge of the column Is determined. The edge of the column 
is oh the vertical line segment extending through the pbcel having this maximal derivative. ^ 

It should be noted that it is not necessary to practice the present Invention in its broadest ssnse that the 

26 left and right edges of each column be found independent of each other. For Instance, as an alternative, one 
of those edges may be found; and then the other edge of the column may be considered as being on the vertical 
line segment parallel to and spaced a preset distance from that found edge, either to the left or to the right 
thereof depending on whether the initially found edge is on the right or left edge, respectively, of the column. 
Once the left and right borders of each column are found, the column detection program then searches 

90 for the bottom border of the column. This search also utilized the fact that, due to the diffraction of light by 
that edge, that edge appears slightly darker on the pixel array than the immediately adjacent areas. I^ore spe- 
cifically, once the x-coordlnates of the two side borders are found, the centerllne of the column is determined 
and used as a reference to locate the third, bottom rectangular area shown In Figure 14. This area, In turn, Is 
used to locate the bottom edge of the column by means of a procedure analogous to the procedures used to 

35 detect the left and right edges of the column. In particular, the gray values in the block are projected onto a 
one-dimenslonal vertical vector. At each point on the vector, a projected value is obtained by adding the gray 
values of all the pixels In the corresponding horizontal line through that point These values are then derived 
and the maximal derivath/e Is obtained. This maximal derivative corresponds to the maximal variation of the 
gray scale, and the bottom edge of the column is considered to be on the horizontal line segment extending 

40 through the pbcel having this maximal variation. 

Once the center line and the bottom point of each column are determined, a smaller window is superim- 
posed over each column. This window which is shown in Figure 16, has the same width as a column; and by 
matching the center and the bottom of the window with those of the column detected above, the window is 
fitted onto the column. Thus, the location of each column Is fully detenmined. 

45 After the column detection program is completed, the feature calculation program begins. As an initial step 
in this program, vark>us reference values are determined for subsequent use; and preferably, a respective set 
of such reference values is determined for each of the two columns illuminated on the pbcel array. More par- 
ticularty, the program selects three reference areas or regions for each column. As shown in Figure 15, one 
region Is located In the bead area inside the column, a second region is located above that bead area, and the 

so third region is located outside the column. 

After the desired reference areas are defined, the program detenmlnes values for the variable Rmodus. Rmax. 
Riniii> Routw* and R«bava- In particular, Rmodw i& s^t equal to the most frequent gray value in the reference region 
In the bead area Inside the column. This reference value is a characteristic of the glass beads in column and 
is used as a threshold value In various subsequent processing steps. In addition, Rm^ Is set equal to the max- 
55 Imum gray value in the reference region in the bead area inside the column, and Rmtn >s ^dt equal to the minimum 
gray value In that inside the column. Router ^ set equal to the average gray value in the reference area outside 
the column, and Rtbcw« Is set equal to the average gray value In the reference area above the glass beads. 
Next, the program begins to extract features related to the reaction that occurred in column 82. The fea- 
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10 



15 



tures extracted include (1 ) parameters related to the cell peltetehape; (2) the red cell agglutinates in the column; 
and (3) the side to side balance of red cells In the column. 

The oeO pellet In the bottom of the column is first obtained by applying a global threshold in the V shape 
region of the column. The default threshold value is 54% of R„edui. In particular, the number of pixels in the V- 
shaped region of the column that have S values less than 54% of R^wiia is determined. The size of the pellet 
Is calculated on the basis of the number of pUels Inside the pellet area, tf the size Is larger than a given number, 
such as 200 pixels, the ceil pellet Is considered significant and the parameters related to the location of the 
cell pellet are determined. These parameters are Illustrated In Figure 16 and are defined as follows: 



Name 


Description 


leftx, lefty 


left-most point in the upper bonder 


rightx, righty 


right-most point in the upper border 


Yllmit 


lowest y-coordlnate In the upper border 


pellety 


bottom-most point In the bottom border 


pelietx 


gravity center of the part below Yllmit 



20 



26 



On the basis of these coordinates, the position of the window Is readjusted by matching the center line of 
the window with pelietx. 

To analyze the pellet shape, the upper border of the cell pellet is fitted with a linear line Y«a+bX This line 
Is determined as follows: 

Assume the upper border of the pellet is defined by a set of points 

(x,.y,).i = 1,N 

where N » rightx " leftx - 3 (The fitting does not includes the four points near the edge on each side). The 
enror of the approximation line is defined as 



N 



35 



40 



To minimize this error, the coefficient a and b have to be 

NSxx - Sx 
a = iSy - bS, 

where 



45 



N 



N 



so 



N 



N 
isl 



The residual of a pellet Is calculated as 
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Residual = 
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20 



26 



30 



35 



40 



40 



SO 



On the basis of the above computation, three feature variable are obtained Including the size of the ceil 
pellet, the eiope of the ceU pellet, and the residual value. These variables are subsequently used to classify 
the agglutinate pattern. 

After these variables are obtained, a fixed mask, shown In Figure 1 7. Is used to cover the whole column 
area. This masic has the same width and shape as a column, and ihe mask is stored in a file and loaded Into 
one memory frame on the Euresys board during Image analysis. By matching the center and the bottom of 
the maak with those of the detected column, the mask is fitted onto the column. 

The next step of the program Is to extract the number of red cell agglutinates and their distribution In the 
column. For this purpose, the bead column is- divided Into five zones, shown in Figure 18, refenred to as the 
positive zone, the negative zone, and Intermediate zones 1, 2, and 3. Generally, the positive zone Is defined 
so as to contain the surface area on the lop of the glass beads, and for example. It may be defined as the area 
above the line Ytop = 80. The negative area Is defined as the cell psllet area in the bottom of the column. If 
there is no cell peUet, then there is no negath/e area. The bead area between the positive and negative zones 
Is divided Into three areas of equal height to form the Intermediate zones 1, 2, and 3. The size of the three 
Intermediate zones Is determined by and Yii„u. and the height. H. of each zone is given by the equation: 



If there Is no cell pellet and, thus, no negathre area, then Yundt defined as 40 pixels above the bottom of the 



The next step In the progriam is to detennine the number of pixels in the positive zone that are illuminated 
at an intensity below a given value, and for example, that given value may be 54% of Rmodw As discussed in 
greater detail below, the number of such pixels in the positive zone is used to determine If a strong poslth/e 
reaction, such as a H reaction, occunned In the column. 

Then, the number of red cell agglutinates located in each of the intermediate zones Is detemilned by means 
of an operation referred to as a top-hat operation and that finds the agglutinates on the basis of the local va- 
riation of the gray values. A top-hat operation is based on two basic operators: dilatation and erosion. The di- 
latation increases the size of an object and the erosion reduces the stee of an object, and an erosion followed 
by a dilatation on an object is usually called as an opening operation. The top-hat transformation of an image, 
denoted h, la defined as: 



where f is the Input image and b is the structuring element function used by the opening (o). In the present 
appilcatfon, the structuring elentent is kernel 7x7, and Figure 1 9 shows an example of the top-hat transforma- 
tion to a line slice of an image. After the dllatlon-eroslon operation, a global cut-off Is applied with a value equal 
to 12% of R„»du»- The number of pbcels in each intermediate zone having gray values, after the top-hat trans- 
formation, greater than 12% of R„»dii» ^ determined. The parameters are then calculated by counting the num- 
ber of pbcels above the cutoff value in the zones 1, 2, and 3. and they oorraspond to the amount of red cell 
agglutinates In these zones. 

The feature calculation program then examines the balance of agglutinates between the left and right 
halves of the colunrm, and in particular, between lower portions of the left and right halves of the column. The 
preferred area of the column that is used to determine that balance is shown in Figure 20; and with reference 
thereto, that area extends upward for a given number of pixels, such as 120 pixels, from the lowest point (Y umft) 
of the upper border of the ceil pellet. The center line of the column is used to separate the area Into two parts, 
as shown in Figure 20. The number and location of red cells In the column and In the cell pellet were previously 
determined during the top-hat procedure and the global threshold, and this data are used to detemiine the num- 
ber of red blood cells on the left and right sides of the column. The balance of agglutinate between the left 
and right halves cf the column is considered as the difference between the numbers of red blood ceils on the 
left and right sides of the column. 

As discussed above, preferably two images, refen'ed to as the front and back images, of each column are 
produced on the pbcel array. The front image is produced, the column is then rotated ISO**, and then the back 
image Is produced. Preferably, values for each of the above-described parameters are obtained for each of 
the front and back images of the column, and then the two values for each parameter are summed. 
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Thus, thefeature extraction program calculatefi thefotlowing features foreach coiunnn: (1) the agglutinated 
red cells In the positive zone; (2) the cell agglutinates In intermediate zone 1, zone 2, and zone 3; (3) the size, 
slope, and residual of the cell pellet and (4) the balance of the agglutinates between the left and right sides. 
These include a total of eight parameters, 

ff Once values for the above-discussed parameters are obtained, the reaction grading program then uses 
these parameter to classify the reaction that occurred in the column into diffenant classes. Generally, the pro- 
gram classifies the reaction as positive or negative; and If positive, as a class +1, +2, +3, or +4 reaction. In 
addition, the program is also capable of identifying a column as having an intermediate reaction, of identifying 
an empty cassette, and of indicating if a cassette cannot be read or if a column cannot be found. 

10 The classifier is a linear decision tree on the basis of the iVIahelanobis distance. This method Is described 
in detail in the boolc "i^ethods Statlstiques da Reconnaissance des fotmosT by O. Qaillat (Ecole Natlonale 
Superieure de Techniques Avancees). Asimpie mathematical definition of different calculations for the meas- 
urement of the separability between classes Is given below. 

15 Discriminate Analysis: . . 

Consider a set of N elements belonging to K classes. Each of the K subsets has Ni, N2, Nk points, and 

is noted as 

20 with n=[1 ..N J and lc=[1..K] and element of subset (I. number of features). 
The center of gravity of this subset Is defined as 



25 




30 and the covariance matrix as 

Ik»2(Xkn-Xk)(X|cn-Xkr 
n=1 

The center of sravity X and the covariance nnatrix T of the whole set are related to those of the subsets 
« ^andTk 



K 



» K K 

T»£Tk*ENk(Xk-X)(Xk-2^) 
k=1 k=1 

56 

Now we define 
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K 

k=1 

which Is the sum of Covarlance matrixes of different classes, termed as the intraciass covariance matrix* and 

K 

W=Tk+2;Nk(Xk-X)(Xk-X/ 



which is the covariance matrix of a set constituted by K points with coefficients Nk, termed as the interdass 
covariance matrix. 

Accordingly we define two additional terms the intraciass variance, v, and the interclass variance, w. as 
20 follows: 

v = trace(V)=2 Z IXkn'^kl^ 
^ k=lns1 

which measures how the elements belonging to one class are grouped around their center of gravity. The small- 
er It is, the nearer the elements are to the center of gravity. If the Intraciass variance was equal to 0, the ele- 
30 mente of a class would be concentrated on their center of gravity. The Interclass variance, w, Is defined as 
follows: 



40 which Is the dispersion between classes. The greater It is, the more classes are separated from each other. If 
ail centers of gravity aggregated to their center X, then the interclass variance would be equal to 0. 

In pattern recognition, the ratio w/V can be used as a measure of separability between classes. If It is great, 
then dassea form compact sets and are all separated from each other. In contrast, if this ratio is small, it will 
be difficult to distinguish classes. 
45 The purpose of discriminate analysis Is to project a set of data Into a J-dimensional sub-space In a way 
that the different classes are most separated. This corresponds to maximize the ratio of the interclass variance 
w over the interclass variance v to be maxinrtum in the sub*space. 

If the data set is projected Into a one-dlmenslonal space with a unitary vector u, then the Intraciass and 
interclass variances of the projected set can be calculated 
50 Vu = uVu 

Wu = u'Wu 

where V Is the Interdass covariance nnatrbc and W Is the intradass covariance matrix. Consider the ratio 

The maximum of this ratio can be found when the gradient 

grad (ru) » ;^x(Vugrad(wJ - WygrBd(vJ) 
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becomes zero. That is, the vector grad(Vu) and grad(Wu) are co-linear. 

grad(Wu) s Xsrad(Vu) 

or 

Wu ° XVu 

5 where X is a parameter. If V can be Inverted, then 

V-m = Xu 

This equation shows that the ratio Tu Is maximum when the direction of the projection axis is defined by the 
eigenvector associated lo the highest eigenvalue of matrix V-^ W. This can also be demonstrated to be true 
for the J-dimensional eub-epace. 
10 To separate the data Into various classes, the feature data is first projected Into a sub-dimensioned space 
to maximize the separability of the data. A classifier has been developed to separate the data in the sub-di- 
mensional space. 

Mahalanobis Distance: 

15 

If a data set consists of K classes (wi, wj,... wj and has Gaussian distribution, then for each point X, the 
probability of X belonging to class Wk can be written as 

fwk(X) = p(WK)-f(X|wO 

with fwkPO is a probability of X belonging to W^; p(W|J is the probability of w^ in .the whole data set; f(X | is 
20 a conditional probability of Xgh^en that It belongs to wic. The classical Bayeslan approachforpattern recognition 
Is to select a class that maximizes 

When the samplea obey the Gaussian distribution, the above approach is equivalent to search the maxi- 
mum of 

gk(X) = - (X - XOTr^pC - XO + Dog(det(T,J) - 21og(p(Wk))] 
25 The first term can t>e interpreted as a square of a distance between X and the center X^, called a Maha- 
lanobis distance. The second term is a correction term dependent on class k but not on X. A set of surfaces 
defined by 

(X - X|JTk-^(X - Xk) = constant 
constitutes a group of concentric eliipaoida.'with center X^, 
so In" the case that the projection space is two-dimensional, the function g^PQ can be written explicitly. The 
calculation is simplified by suppressing the corrective tBtm If the equlprobaWllty of the classes is considered. 

Assume the sample X is represented by the coordinates (x,y), the center of gravity of class k is pCk, Y|J. 
and the inverse of the covarlance matrix of class, kTtr\ Is 



which Is symnwtrk! Tit.i2'=Tk^t. Thus, assume Gk(x.y) « -gk(X) and It can be written as 
Gk(x,y) » Tk,ii(x - Xk)2 + 2xTk.i2(x - Xk)(y - yic) ♦ Tk^(y - yk)^ 
45 Gk(x,y) Is a second degree polynomial, and we want to find the minimuhi Gk. 

The samples are graded into different classes on the basis of a linear decision tree. This program classifier 
first separates the samples Into two principal groups. One group, Group I. contains the classes +1, +2, +3, 
and +4 reactions, and the other group, Group II. includes 0 and very weak positwe reactions. The classifier 
identifies these following classes: 
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Class 


Description 


0 


Negative reaction 


1 


Positive reaction 


2 


Positive reaction 


3 


Positive reaction 


4 . 


Strong positive reaction 


-2 


Can't read cassette or column not fbund 


-4 


Intermediate reaction 


-5 


Empty casGfitta 



Figure 21 shows a flow chart summarizing the global structure of the classifier. The feature data used in 
each decision are shown in the Figure and are detailed as follows: 

Empty Column: To determine if a column is empty, the total red cells In a whole column, including the negative 
20 zone, positive zone, and three intermediate zones, are summarized. If the value Is less than a given number 
such as 500, the column Is graded as an empty column. That is 

Sum^whole = PPos + PNeg + Zonal + Z6ne2 ♦ Zone3 
If (Sum_whole<500), then {empty column} 
Separation into two groups: The negative class reactions are aspareted from the majority of the positive class 
26 reactiona on the basis of the eize of the cell pellet and the sum of the agglutinates In zones 1 , 2, and 3. That 
Is: 

if (PNeg<50D) or (Sum_Zones>800). then positive reaction {1 . 2. 3. 4} 
else negaUve reaction {0,1} 
Classification of {1,2,3,4): The positive reactions are further separated Into dass +1 , +2, +3, and +4 reactions 
30 on the basts of the distribution of agglutinates in the column. The distribution of the agglutinates Is represented 
by the five features: PPos, Zone1. 2one2, Zone3, and PNeg. Table 1 lists the average feature data of four 
types of positive reactions. The results indicate how the feature data vary among positive reactions. 



Table 1: 



3S 



40 



45 



Average feature data of different reaction classes. 


Features 


Class 1 


Class 2 


Class 3 


Class 4 


PPos 


22 


304 


2183 


3860 


Zone 1 


358 


1238 


1745 


259 


Zone 2 


302 


1034 


700 


43 


Zoned 


359 


786 


270 


58 


PNeg 


2085 


764 


26 


19 


Number of samples 


281 


328 


173 


45 



The classification is done by Mahalanobis distance, vkrtiich is able to separate correctly the different types 
of reactions. 

Classification of (0.1): The separation of these two classes uses the following features: (I) the sum of agglu- 
tinates In zones 1 , 2, and 3, (II) the slope and (ill) the residuals of the ceil pellet upper border, and (Iv) the side 
to side balance of red cells. In order to Increase the separability between the two classes, the discriminate 
analysis la applied. The above four features are projected Into one dimension space. The variation of the sep- 
arability Is significant Before the projection, the separability is equal to 0.70 and after projection, 1.38, This 
means that the elements of each dass are more regrouped around their center of gravity and the distance be- 
tween the two classes Is Increased. The projection Is calculated as following: 
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Pfojection » Q.78xAgglutinatd8 + 0.48xre6idual8 -i- 0.36Xs)op6 4- 0.13xbatance 
The coeff iclenta are calculated by discriminate analy8i& They ahow the Importance of each feature for 
the projection. 

The limit k)6tween claaaea is estimated by Mahalanobia distence. The result may also include an indeter- 
5 minate class. It may be defined by the difference between the two Mahatanobis diatences calculated for both 
classes, if the absolute difference Is below 1.2. then the sample Is classified as a questionable reaction. This 
third class Is on the boundary between classes 0 and 1 . 

While It Is apparent that the Invention herein disclosed Is well calculated to fulfill the objects previously 
stated. It will be appreciated that numerous modifications and embodlmente may be devised by those sicllled 
io in the art, and it is intended that the appended claims cover all such modifications and embodiments as fall 
within the true spirit and scope of the present invention. 



Claims 

15 

1. A method for analyzing a solution for an agglutination pattern, comprising: 

producing an Illuminated Image of the solution on an array of pixels; 

assigning to each pixel in the array and on said Illuminated Image, a data value representing an 
intensity of the illuminated image on the pixel; and 
20 processing the data values according to a predetermined program to determine If an agglutination 

pattern is present In the solution: and If an agglutination pattern Is present, classifying the agglutination 
pattern into one of a plurality of predefined classes. 

2. A method according to aalm 1, wherein the processing step includes the steps of: 
26 separating the array of pbceis Into a plurality of zones; 

processing the date values for the pbceis in each zone according to a respective predetermined pro- 
cedure to determine values for a predefined set of variables; and 

processing said determined values to determine whether an agglutination pattern is present and, 
If an agglutination pattern Is present, to classify the pattern Into one of the predefined classes. 

^ 3. A method according to Claim 2, wherein the solution Is in a column having a bottom, a lower portion, and 
a filter disposed in said lower portion, and wherein the separating step includes the steps of. 
identifying a first zone of pixels adjacent a top of the filten and 
Identifying a second zone of pUels adjacent the bottom of the column. 

^ 4. A method according to Claim 2 or Claim 3, wherein the step of processing the date values to determine 
values.for the predefined set of variables Includes the steps of: 

assigning to a first variable a value reprssentlng the number of pixels In the first zone having a 
date value less than a first value; 

assigning to a second variable a value representing the number of pbceis In the second zone havi ng 
a date value less than a second value. 

5. A method according to Claim 4, wherein 

the step of processing the date values to detemilne values for the predefined set of variables further 
Includes the step of esteblishlng a reference value; 
^ the step of essigning a value to the first variable includes the step of assigning to the first variable 

a value representing the number of pbceis In thef Irst zone having a date value less than a given percentege 
of the reference value; and 

the step of assigning a value to the second variable includes the step of assigning to the second 
variable a value representing the number of pbceis in the second zone having a data value less than a 
^ gh/en percentage (}f the reference value. 

8. A method according to Claim 5, wherein the establishing stop Includes the steps of: 
defining a reference area on the pbcel anray; 

establishing as the reference value, the most common data value for the pbceis In the reference 

5ff area. 

7. A method according to Claim 6. wherein the defining step Includes the step of defining the reference area 
In the Image of the column on the pbcel array. 
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8. A method according to Claim 7, wherein the step of defining the reference area in the image of the column 
further Includes the step of defining the reference area in the Image of the filter on the pixel array. 

9. A method according to Claim 4, wherein the separating step further includes the step of Identifying a plur- 
ality of additional zones of pixels between said first and second zones. 

10. A method according to Claim 9, wherein the step of processing the data values to detarmine values for 
the predefined set of variables further includes the step of assigning to an additional set of variables values 
representing the numbers of pixels in the additional zones having data values less than a third value. 

11. Amethod according to anyone of Claims 1 to 10 whersin the assigning step includes the step of identifying 
the pixels on the array and in the illuminated image of the solution. 

12. A method according to Claim 11, wherein the illuminated Image has left and right edges, and the identi- 
fying step includes the step of identifying the pUels on the left and right edges of the illuminated image. 

13. A method according to any one of Claims 1 to 12 wherein the solutton is In a column having a bottom, 
and wherein: 

the step of processing the data values for the pixels include the step of proceasing the data values 
to determine if a pellet of cells is located in the bottom of the column; and 

the predefined set of variables Include the size of the pellet; and the stope and shape of a top sur- 
face of the pellet 

14. Amethod according to aalm 1 3, wherein the column Includes left and right halves, and each of said halves 
includes a number of btood cells, and wherein the predefined set of variables further Includes the differ- 
ence between the numbers of blood cells on the left and right habeas of the column. 

15. A method for analyzing a solution for an agglutination pattern, comprising: 

producing an Illuminated image of the aoiutlon on an array of pixels; 

assigning to each pixel on the Illuminated Image, a data value representing the intensity of the il- 
luminated image on the pixel; 

separating the array of pixels into a plurality of zones; and 

processing the date values of the pixels in each region according to a predetermined program to 
detonnlne whether an agglutination pattern is present In the solution. 

16. Amethod according to Claim 15, further including the step of, if an agglutination pattern is present in the 
solution, classifying the agglutination pattern Into one of a plurality of predefined classes. 

17. A method according to Claim 16, wherein the processing step includes the steps of: 

processing the date values for the pbcels in each zone according to a respective predetermined pro- 
cedure to determine values for a predefined set of variables; and 

processing said determined values to determine whether an agglutination pattern is present 

1 6. A method according to any one of Claims 1 to 1 7 wherein the solution is in a contdner leaving a bar code 
Identifying date, and further comprising: 

reading the bar code on the cassette; and 

omitting the processing step if the bar code does not identify all of preselected date. 

19. A method according to any one of Claims 1 to 18 wherein the solutton is in a container, and further com- 
prising: 

if the container has preselected alignment marks, positioning the container to locate the alignment 
marks in preset positions to align the conteiner with the pbcel array, and omitting the processing step if 
the alignment marks are not located in the preset positions; and 

if the conteiner does not have preselected alignment mailcs, omitting the processing step. 

20. A system for analyzing a solution for an agglutination pattern, comprising: 

nneans for holding the solution; 
a pixel array; 

illumination means for producing an Illuminated Image of the solution on the pb(6l an-ay; 

means to assign to each pbcel on the Illuminated Image, a date value representing the Intensity of 
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the illuminated Image on the pixel; and 

processing means to process said data values according to a predetermined program to determine 
whether an agglutination pattern is present In the solutton In the sample contalnen and If an agglutination 
pattern is present, to classify the agglutination pattern Into one of a plurality of predefined classes. 

6 

21. Asystem according to Claim 20, wherein the holding means includes 

a base; 

a frame fbr holding the solution and pivotaliy connected to the t>ase; 

moving means connected to the frame to pivot the frame and the solution to position first and sec- 
' end opposite sides of the solution facing the pixel array. 

22. A system according to Claim 21 , wherein the solution is contained within a cassette, and the frame forms 
an elongated channel for holding the cassette. 
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23. Asystem according to Claim 22, wherein the elongated channel Inciudes first and second opposite, open 
ends to allow sliding movement of the cassette into, through and from the elongated channel. 



24. A system according to any one of Claims 20 to 23, further comprising 

a storage rack located adjacent the frame, and fonnning a multitude of elongated channels for hold- 
ing a multitude cassettes; and 
^ Indexing means connected to the storage rack to move the storage rack through a series of posi- 

tiona to align each of the elongated channels of the storage rack with the elongated channel of the frame. 

25. Asystem according to any one Claims 20 to 24, wherein the means for holding the solution Inciudes bar 
codes Wentlf ying the holding means, and the system further Includes a code reader for reading the bar 

26 code and identifying the holding means. 

26. A system according to any one of Claims 20 to 25, wherein the mesns for holding the solution include at 
least one positioning mark, and further comprising means to locate the holding means In the frame with 
the positioning marks in predetermined locations. 

so 

27. Asystem according to any one of Claims 20 to 26, the processing means includes: 

means to separate the array of pixels into a plurality of zones: 

means to process the data values for the pixels In each zone according to a respective predeter* 
mined procedure to determine values for a predefined set of variables; and 
^ means to process said determined values to detemiine whether an agglutination pattern is present, 

and if an agglutination pattern Is present to classify the pattern into one of the predefined classes. 

28. A system according to Claim 27, wherein the solution is in a column having a bottom, a lower portion, and 
a filler disposed In said lower portion, and wherein the means to separate the array of pixels includes 
means to Identify a first zone of pixels adjacent a lop of the filter, and to Identify a second zone of pixels 

■ adjacent the bottom of the column. 



29. Asystem according to Claim 28. wherein the processing means further Includes 
means to define a reference area on the pixel array; and 

means to establish as a reference value, the most common data valuator the pixels in the reference 

area. 

30. Asystem according to any one of Claims 20 to 29, wherein the means to assign a data value to each pbcel 
includes means to identify the pbceis on the anray and In the illuminated image of the solution. 

w 31 . A method for analyzing a solution for an agglutination pattern, comprising: 
producing an illuminated Image of the solution oh an amy of pbceis; 

assigning to each pUel on the illuminated image, a data value representing the intensity of the il- 
luminated image on the pbcel; 

processing the data values according to a predetermined program to determine whether an agglu- 
w tinatlon pattern is present In the solution and, If an agglutination pattern is present, to Uentify the pattern 

as belonging to one of a multitude of predetennnlned classes, the processing step Including the steps of 
1) processing the date values to determine If cells are present in the column, 
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ii) processing the data values according to a first routine to classify the solution preliminarily as being 
In one of two groups, 

ill) If the solution Is classified In a first of the two groups, then processing the data values according 
to a second routine to determine whether an agglutlnatton pattern is present In the solution, and 

if the solution is classified in a second of the two groups, then processing the data values according 
to a third routine to classify the agglutination pattern Into one of a first set of the predefined classes. 

32. A method according to Claim 31 , wherein the predetermined classes Include: a negative reaction, class, 
and a weak poslth^e reaction class, and the step of processing the data values according to the second 
routine includes the step of distinguishing between the solution belonging In the negatWe reaction class 
and the weak positive reactton class. 

33. Amethod according to Claim 32, wherein the predetermined classes further include other positive reaction 
dasses, and the step of processing the data values according to the third routine Includes the step of dis- 
tinguishing between the solution belonging In the strong positive reaction class and the other positive re- 
action classes. 
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Image Acquisition 

(1) Initialize Camera and Euresys Board; 
(Z) Set Euresys Board Parameters; 

(3) Grabe Images; 

(4) Read &Wrlte images. 



Column Detection 

( 1 ) Determine the Location of the Columns in 

an Image; 

(2) Generate window to cover column 



Feature Calculation 

( 1 ) Extract Red Celt Agglutinates; 

(Z) Extract Ceil Pellet and Analyze theQQDQQ 

Pellet Shape; 
(3) Estimate the Balance of Red Cells in a 
Column, 



Reaction Grading 

( 1 ) Combine the Features; 

(2) Grading the Reaction on the Basis of 
the Feature Combination* 
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